Quantum-classical correspondence in the wavefunctions of Andreev billiards 
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We present a classical and quantum mechanical study of an Andreev billiard with a chaotic normal 
dot. We demonstrate that in general the classical dynamics of these normal-superconductor hybrid 
systems is mixed, thereby indicating the limitations of a widely used retracing approximation. We 
show that the mixed classical dynamics gives rise to a wealth of wavefunction phenomena, including 
periodic orbit scarring and localization of the wavefunction onto other classical phase space objects 
such as intermittent regions and quantized tori. 

PACS numbers: 74.45.+C, 75.45.+j, 03.65.Sq 



Wavefunction phenomena in closed and open quan- 
tum dot systems have attracted much attention in recent 
years. Besides the ongoing theoretical interest in a simple 
yet detailed description of scarring0,|2], a range of phe- 
nomena have been studied including the connection be- 
tween conductance oscillations, transmission resonances 
of ballistic semiconductor dots and microwave cavities, 
light emission from dielectric cavities and the localization 
of the wavefunction onto classical phase space objectsQ 
such as Kolmogorov-Arnold-Moser (KAM) isiands£|, hi- 
erarchical regions and periodic orbits It has been 
also pointed out, that certain states of a closed quan- 
tum dot, associated with KAM islands or bouncing-ball 
trajectories, actually survive the coupling of the dot to 
external leads, since level broadenings are not uniform, 
and this can give rise to measurable transport effectsQ. 

Systems consisting of a ballistic quantum dot cou- 
pled to a superconductor, which are commonly called 
Andreev-billiards (ABs)ja, 0] raise new questions of 
quantum-classical correspondence addressed by only a 
few works |H [H E H beforehand. The key phys- 
ical process taking place in normal-superconductor (N- 
S) hybrid nanostructures is the Andreev reflection [l^J, 
whereby electron-like quasi-particles with energies Ej? -(- s 
(Ep is the Fermi energy) are coherently scattered into 
Fermi sea holes of energy E-p — e (and vice versa) at the 
N-S interface if e is smaller than the superconducting 
gap A. For e = the Andreev reflected electrons (holes) 
perfectly retrace their classical trajectories as holes (elec- 
trons). For e > the velocities of the quasiparticles are 
not exactly reversed because of the 2e difference in their 
energies, but this is not considered in most of the litera- 
ture. 

The validity of this widely used retracing approxima- 
tion is questionable for ABs with chaotic normal dots, if 
the time between two Andreev-refiections is significantly 
longer than the ergodic time, because the small misalign- 
ment between the velocities of the electrons and holes at 
the N-S interface leads eventually to divergent trajecto- 
ries. One notable exception, in which the consequences 
of the non-exact velocity reversal were investigated is the 



paper by Silvestrov et al. yj|. They studied an Andreev- 
billiard with a chaotic normal dot and found that even 
though the motion in the normal dot could be character- 
ized by non-zero Lyapunov-exponents, the existence of 
an adiabatic invariant in the N-S system to good approx- 
imation confines the electron-hole orbits to tori. Never- 
theless, the exponential divergence of nearby trajectories 
does manifest itself through a gap in the low energy den- 
sity of states (DOS) and through the eigenfunctions of 
the system, which were predicted to exhibit a peculiar 
localization property, different from those known in nor- 
mal billiards. However, no quantum calculations were 
performed to support the latter finding. 

Recently, Wiersig 0] has shown that the diffraction 
occurring at the points separating the normal and super- 
conducting segments of the dot boundary (coined "criti- 
cal points") also plays an important role in the classical 
dynamics of ABs. Electrons (holes) hitting the boundary 
at nearby points belonging to different types of segment 
will undergo cither Andreev or normal reflection and will 
be scattered into different quasi-particle states (ie into 
holes or electrons). Owing to the special geometry this 
diffraction effect was not observed in Ref. 

In this Letter we show, that in a more general case, 
the interplay of the non-exact velocity reversal and the 
diffraction at the critical points leads to the breakup of 
the (adiabatic) tori in at least certain parts of the phase 
space [see Fig.^a),(b)] rendering the classical dynamics 
mixed. While it has been known for a long time, that 
an applied magnetic field can render the ABs (weakly) 
chaotic d, Q, to our knowledge the possibility that the 
dynamics of ABs can be irregular even for zero magnetic 
field has not yet been addressed in the literature. We 
argue moreover that the combination of these two effects 
results in a wealth of phenomena not only in the clas- 
sical but also in the quantum dynamics of ABs. What 
makes these N-S hybrid systems additionally interesting 
and novel in the realm of quantum chaos studies is that 
the wavefunction is spinor and not scalar, as in normal 
billiards. We calculate, for the first time, the exact quan- 
tum eigenfunctions of a two dimensional AB with chaotic 
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normal dot. Using the Wigner transform of the eigen- 
states we find that most of them can be associated cither 
with the regular or the chaotic regions of the energy sur- 
face. 

The particular system used in our numerical calcu- 
lations is the Sinai- Andreev (SA) billiard^]- It con- 
sists of a Sinai-billiard shaped normal dot and an at- 
tached (infinite) superconducting lead (see Fig. ^c) for 
the geometry). First, we consider the classical dynamics. 
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FIG. 1: (color online) The Poincare section of the SA billiard 
for h = 0.58 (a) and h = (b). (For other parameters see 
[22]]). Each dot represents a starting point of an electron tra- 
jectory; the red and blue dots correspond to very elongated 
regular islands. The torus shown with red dots in (a) is pro- 
jected onto real space (c). 

Fig. ^a,b) shows the Poincare section (PS) which is de- 
fined as follows. The coordinates of each point in PS are 
the position y of the electron starting from the N-S in- 
terface and the tangential velocity component v y in units 
of v e = vp^Jl + e/ E-p, where vf is the Fermi velocity (y 
is measured from the lower corner of the N-S interface). 
As can be observed, the PS resembles the Poincare maps 
of generic normal systems with mixed classical dynam- 
ics. In large regions of the energy surface the motion 
is chaotic, while islands of stability are preserved only 
around such unstable periodic orbits of the isolated nor- 
mal system, which hit the N-S interface at right angle. 
We emphasize that in contrast to Ref. these phase 
space structures are regular islands not only in the adia- 
batic approximation. Depending on the instability of the 
orbit and the value of e/Ep. the presence of the supercon- 
ductor can indeed stabilize the motion, as can be checked 
by computing the stability matrix of the corresponding 
electron- hole orbit. Most of these islands are extremely 
elongated (see eg Fig.^a,b), where they are denoted by 
red and blue dots). Moreover, considering the bunch of 
electron trajectories lying on a particular torus, its pro- 
jection onto the real space is "squeezed down" away from 
the N-S interface [Fig. H[c)] . Both phenomena are con- 
sequences of the exponential divergence of nearby trajec- 
tories in the nor mal dot [Toll. For islands centered on the 
least unstable orbits (eg the one at h = 0) these effects 



are less pronounced [see the large island in Fig. ^b)] . 

Furthermore, strips of regular, intermittent-like mo- 
tions can also be observed in Fig. ^a), in particular 
around v y /v c ps 0.7, 0.3, —0.2, and —0.5. These re- 
gions correspond to such initial conditions for which the 
electron and hole trajectories will only be slightly dif- 
ferent, but this nearly periodic orbit will slowly drift 
in phase space0, |5fl| [see also Fig. Hfb)]. This dy- 
namics resembles the intermittent behavior in normal 
billiards|23j where this near-integrable motion usually 
evolves in the vicinity of tori or isolated periodic orbits. 
In the present case [Fig. 12a)], however, it is a conse- 
quence of the retroreflection mechanism and is therefore 
a peculiarity of ABs. Nevertheless, owing to this drift 
either the electron or the hole trajectory will eventually 
reach the edge of the N-S interface and then by hitting 
the normal wall instead of the superconductor, escapes 
to the chaotic sea. 

In what follows, we show how the properties of the 
classical phase space leave their fingerprints on the eigen- 
states of Andreev billiards. First, we briefly summarize 
the quantum treatment of the system. The spinor wave- 
function ip(r) — [u(r), v(r)] T (with u(r) electron and v(r) 
hole components) of the N-S hybrid systems satisfies the 
Bogoliubov-de Gennes equations T-Lip^r) = eip(r), where 
H = n Q a z + A(r)a x , and H = -h 2 /2mV 2 - is the 
single-particle Hamiltonian with Dirichlct boundary con- 
ditions at the normal walls. Here \i is the chemical poten- 
tial and o~ z , o~ x are Pauli matrices. We assume that the 
superconducting pair potential is A(r) = A constant 
inside the lead and zero in the N region^. The calcu- 
lation of the wavefunction requires two steps: first we 
obtain the energy levels of the SA billiard using a quan- 
tum mechanically exact secular equation which can be 
derived invoking the scattering approach of Ref. Q . The 
method also furnishes us with the wavefunction in the 
superconducting lead and at the N-S interface. Then the 
boundary integral method is employed to find the wave- 
function in the normal dot|l8l|. We work in the regime 
Sn -C Et <C A <C -EfU where 8 n is the mean level 
spacing of the isolated normal dot and Et is the Thou- 
less energy ji^ ■ 

The classical-quantal correspondence of phase space 
structures and eigenstates of a given system can be stud- 
ied with the help of the Wigner function[2l]]. In order 
to compare the Wigner function with the PS, we calcu- 
lated the projection of the Wigner function[2j| of both 
the electron and the hole components of the wavefunction 
onto the PS: 

Wp(y,P,)-^ JdY e -^ Y C{x,y-^)ax,y+\) (1) 

evaluated at the N-S interface (ie x = 0), ( is either 
u(r) or v(r) and p y is the parallel (to the interface) com- 
ponent of the momentum. The Wigner function is not 
positive definite and usually exhibits rapid oscillations 



that can obscure the physical content. For this reason, 
like in Ref. j2(|, we smoothed the projection Wp with a 
Gaussian, which was chosen narrower than the minimum 
uncertainty Gaussian. Note that Wp is symmetric in p y , 
since the system is time reversal invariant. This sym- 
metry is absent in the Poincare map, because incidental 
electrons (holes) are not taken into account. 

In common with normal billiards, most of the eigen- 
states of the SA billiards can be classified as 'chaotic ' or 
'regular' [3 El- The chaotic eigenstates can further 
be subdivided into two groups: i) crgodic-like, and ii) 
scarred states. Regarding the regular states, they either 
can be associated with the intermittent-like regions of 
the phase space or with quantized tori. We now consider 
each case. 

For ergodic-like eigenstates (an example is shown in 
Fig. Ufa)) both the electron and the hole components 




FIG. 2: (color online) The probability densities of the elec- 
tron (|m| 2 ) and hole (\v\ 2 ) components of a chaotic eigenstate 
(e/Ao = 0.060395) of the SA billiard (a). The corresponding 
smoothed projections of the Wigner function onto the PS (b) 
in units of y/W and p y /p e , Py/ph for the electron and hole 
components respectively. Herep e = mv e andph(e) = p e (— s). 
Eight equally spaced positive contours (red lines) and five 
negative contours (green lines) are plotted. The distribution 
of the scaled probabilities z e (red line), Zh (green line) and 
the Porter-Thomas distribution (blue line) (c). 

of the wavefunction seems to cover the normal dot in 
a roughly uniform way. Nevertheless, they display dif- 
ferent interference pattern which translates also into the 
corresponding projections of the Wigner functions. Ex- 
amining Fig. Ufa) and Fig. Efb) one can see that Wp 
has high amplitude in such regions, which corresponds to 
chaotic regions in the PS. This observation is reinforced 
by the good agreement between the probability distri- 
bution of the scaled local densities z e (r), Zh(r) and the 
Porter-Thomas distribution V(z) = l/v / 27rzcxp(— z/2) 
shown in Fig. |2Jc). Here z e (r) = (A/N e )\u{v)\ 2 , where 
Af e = J |«(r)| 2 d 2 r, and the integration is performed over 
the area A of the normal dot (an analogous definition ap- 
plies for the hole density z/,(r)). Both the electron and 
the hole components of chaotic eigenstates of ABs can 
thus be considered, similarly to the eigenstates of chaotic 
normal billiards |2 ll l22j . as being a superposition of infi- 
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nitc number of plane waves with fixed wave number, but 
with random directions and amplitudes. 

The old wisdom of the quantum chaos: 'the scars are 
scarce' seems to hold also for ABs. Out of 99 eigenstates 



b) 




FIG. 3: (color online) The probability densities of the elec- 
tron (|m| 2 ) and hole (\v\ 2 ) components of a scarred eigenstate 
(e/Ao = 0.88323) of the SA billiard (a) and the correspond- 
ing unstable periodic orbit (b). The probability densities of 
a decoupled state at e/Ao = 0.42901 (c) and the marginally 
stable periodic orbit family over which the hole component 
shows enhancement (d). 

of the system with geometrical parameters listed in the 
caption of Fig. ^a), only 2 can be classified as scarred, 
one of them is shown in Fig.^a) along with the periodic 
orbit that scars the electron component (Fig. Efb)). We 
call these states 'genuine Andreev-scarred' to distinguish 
them from the other type of states which are decoupled 
from the superconductor [l^ and in certain cases also 
show enhancement over periodic orbits (see Figs. |3{c) 
and (d)). In the case of decoupled eigenstates the am- 
plitude of the wavefunction at the N-S interface is small 
and the probability of finding the quasiparticlc in cither 
the electron or hole state inside the dot is enhanced. For 
the example shown in Fig. Old), integrating |w(r)| 2 over 
the area of the billiard one finds that this probability is 
91% whereas only 7% for the electron component (and 
2% corresponds to the probability of quasiparticles in the 
superconductor). In contrast, for the eigenstate shown in 
Fig- 01 a ) the wavefunction has an apparently finite value 
at the N-S interface and the probabilities arc 56% (elec- 
trons) and 31% (holes). 

Regarding the regions of regular motion, we found that 
both the intermittent-like parts of the phase space and 
the regular islands can give rise to quantum eigenstates. 
The electron and the hole components of regular eigen- 
states display very similar interference patterns and the 
probabilities of the electron and hole state in the dot have 
close values. Furthermore, the localization of both the 
eigenfunction and the corresponding smoothed Wp onto 
the underlying phase space structure can be observed. 
As for the islands of stability shown in Fig. Qfa), they 
enclose a very small area and thus for the numerically ac- 
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ccssiblc wavelengths they are not resolved quantum me- 
chanically. A much more amenable system for studying 
the correspondence between quantum eigenstates and a 
regular island can be obtained by attaching the super- 
conducting lead at h = 0, since in this case the emerging 
regular island is larger [see Fig.^b)] . Our findings on the 





FIG. 4: (color online) The probability density of the electron 
(|m| 2 ) and hole (|t'| 2 ) components of an intermittent eigenstate 
(e/A = 0.44652) of the SA billiard (a). One of the slowly 
drifting electron-hole orbits during the intermittent-like mo- 
tion (b) which form the regular strip around v y /v e ~ —0.2 in 
Fig.^a). The orbit comprises electron (red) and hole (green) 
trajectories. The smoothed projections Wp for the electron 
and the hole components of the eigenstate (c). Eight equally 
spaced positive contours (red lines) and one negative contour 
(green line) are plotted. The blue curve marks the location 
of that region of the PS, in which the drifting electron- hole 
orbits, like that one in Fig. |3Ib), stay longer than r H (see 
text). The electron component of a wavefunction localized 
onto a quantized torus when h — (d). The hole component 
is almost identical and thus not shown here. 

eigenfunction properties listed above are based on study- 
ing this case [see also Fig. Eld)], which is, however, not 
treated here in more detail for the lack of space. Finally, 
we found that as expected, those intermittent-like regions 
support quantum eigenstates, which correspond to stay- 
ing times longer than the Hcisenberg time r H = 2Ti/8^- 
An example is shown in Fig. Ufa), while Fig.^Jc) shows 
that the smoothed Wp also has high amplitude in the 
corresponding region. 

In summary, we investigated the quantal-classical cor- 
respondence in Andreev billiards with a chaotic normal 
dot. We showed that in contrast to the common belief, 
the interplay of critical points and non-exact velocity re- 
versal can render the classical dynamics non-intcgrablc 
even for zero magnetic field. It was also shown that the 
eigenstates of the system can be classified as chaotic or 
regular corresponding to different regions of phase space. 
This implies that while the retracing approximation has 
been proved to be useful in understanding the energy de- 
pendence of the density of states [l3L Il7j. it may not be 
adequate when addressing the properties of individual 
eigenstates. Experimentally, the eigenstates of the An- 
dreev billiards might be studied using scanning tunneling 
probe, and those which comprise quasiparticlc density en- 



hancement in certain part of the dot should be discernible 
from the ergodic ones. 
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